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ABSTRACT 

The interaction of relativistic magnetized ejecta with an ambient medium 
is studied for a range of structures and magnetization of the unshocked ejecta. 
We particularly focus on the effect of the ambient medium on the dynamics of 
an impulsive, high-sigma shell. It is found that for sufficiently high values of 
the initial magnetization cxo the evolution of the system is significantly altered 
by the ambient medium well before the shell reaches its coasting phase. The 
maximum Lorentz factor of the shell is limited to values well below ctq; for a 
shell of initial energy E = erg and size ro = lO^^Tso cm expelled into 

a medium having a uniform density Ui we obtain Fmax — ^80{E52 /T^o'n.i)^^^ in 
the high sigma limit. The reverse shock and any internal shocks that might 
form if the source is fiuctuating are shown to be very weak. The restriction 
on the Lorentz factor is more severe for shells propagating in a stellar wind. 
Intermittent ejection of small sub-shells doesn't seem to help, as the shells merge 
while still highly magnetized. Lower sigma shells start decelerating after reaching 
the coasting phase and spreading away. The properties of the reverse shock then 
depend on the density profiles of the coasting shell and the ambient medium. 
For a self-similar cold shell the reverse shock becomes strong as it propagates 
inwards, and the system eventually approaches the self-similar solution derived 
recently by Nakamura & Shigeyama. 



1. Introduction 



The interaction of relativistic ejecta with the surrounding medium is an issue of con- 
siderable interest. During the early stages of the evolution a double shock structure forms, 
consisting of a forward shock that propagates in the ambient medium, a reverse shock cross- 
ing the ejecta and a contact interface separating the shocked ejecta and the shocked ambient 
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medium. Under certain conditions, an observable flash of electromagnetic radiation is ex- 
pected to be emitted during the propagation of the reverse shock. In the fireball scenario 
commonly adopted, the naive expectation has been that optical flashes associated with re- 
verse shock crossing should be quite common (e.g., Kumar & Panaitescu 2003; but c.f. 
Nakar & Piran 2004), however, despite considerable observational efforts only a few have 
been detected, indicating that such flashes are rare. 

It has been proposed that the paucity of optical flashes may be attributed to an early 
onset of a R-T instability (Levinson, 2010a,b), or strong magnetization of the ejecta (e.g., 
Zhang & Kobayashi 2005; Mimica et al. 2009). The latter is anticipated if the free energy 
is extracted magnetically in the form of a Poynting-flux- dominated flow (e.g., Levinson & 
Eichler 1993; Lyutikov & Blandford 2003 ; Giannios & Spruit 2005). Stationary magnetic 
outflows allow, in general, only partial conversion of magnetic-to-kinetic energy, implying 
high magnetization at the onset of the afterglow phase. A better conversion can be achieved 
if the outflow is coUimated into a small opening ang le, e ~ V-^ (Komissarov et al. 2009), 
though corking from a star, as in the coUapsor model for GRBs, can alleviate the latter 
condition (Tchekhovskoy et al., 2010; Komissarov et al. 2010). However, even then cr ~ 1 is 
anticipated at best (Lyubarski 2009). 

Recently it has been shown (Granot et al., 2010; hereafter GKSIO; Lyutikov 2010a,b) 
that time-dependent effects may play a crucial role in the acceleration of a magnetized 
flow. These authors considered the acceleration of a spherical, impulsive high sigma shell 
of initial width A = ro and magnetization o"o, expelled by a central source. They have 

1/3 2/3 

shown that, unlike a stationary flow for which acceleration ceases at Foo ~ ctq , a^o ~ ctq ' 
the impulsive shell continues accelerating even after loosing causal contact with the central 
source until reaching nearly complete conversion of magnetic energy into bulk kinetic energy. 
The terminal Lorentz factor of the shell is ~ ctq. During the acceleration phase, that 
they term "magnetic rocket acceleration", the major fraction of the shell energy is contained 
in a layer of width 2ro, bounded between the front of a rarefaction wave reflected from the 
central source and the head of the shell. The average Lorentz factor of the shell, roughly 
equals the Lorentz factor of the fluid at the rarefaction front, evolves as < P >oc t^^^. The 
structure of this layer is well described by a self-similar solution. Once the shell enters the 
coasting phase its width starts growing and its magnetization continues to drop. 

In this paper we consider the interaction of the shell with the external medium and show 
that in the high-sigma limit the evolution of the system is dramatically altered. A preliminary 
account of the effect of the ambient medium is given in GKSIO, and a more detailed discussion 
in Lyutikov (2010a,b). Speciflcally, it is shown that for initial magnetization cro larger than 
some critical value, deceleration of the contact interface commences well before the shell has 



- 3 - 



reached the coasting phase, when it is still highly magnetized. The Lorentz factor of the 
contact discontinuity evolves as T^, oc t^^/^, while the rear boundary of the unshocked shell 
is still accelerating. The maximum Lorentz factor of the compressed shell is then limited 
to r << (To- The reverse shock in this case is very weak or nonexistent, and within the 
framework of ideal MHD no internal dissipation is practically expected (by either the reverse 
shock or any internal shocks that might form in a multi-shell scenario) before the onset of 
the afterglow phase. Lower sigma shells start decelerating after reaching the coasting phase. 
The properties of the reverse shock then depend on the density profiles of the coasting shell 
and the ambient medium, as discussed in §2.2[ 




shock 



Fig. 1. — Schematic representation of the double-shock system. There are three characteristic 
surfaces: a forward shock propagating in the ambient medium, a reverse shock sweeping the 
ejecta, and a contact discontinuity separating the shocked ejecta and the shocked ambient 
medium. The Lorentz factors of the three surfaces, measured with respect to the unshocked 
ambient medium, are indicated. Quantities in the shocked ambient medium and shocked 
ejecta are denoted by subscripts 1 and 2, respectively. 



2. Thin Shell Model 

We consider the interaction of a cold, magnetized shell with an ambient medium having 
a density profile pi{r) = air~''. The ejecta is characterized by a Lorentz factor 7e(r, t), 
density Pe{r,t) and magnetic field vector 6^(r, t), assumed to be given, where the 4- vector of 
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the magnetic field, b^, is defined in appendix |X1 To simplify the analysis we shall assume 
a spherical shell with a purely toroidal magnetic field, viz., 6^ = (0,0,0,6e). The magnetic 
pressure of the shell is then 6^/2 and the corresponding sigma parameter is Ue = b^/pe- The 
structure of the shocked shell is shown schematically in Fig. [TJ The subscript 1 refers to 
the shocked ambient medium and 2 to the shocked shell. The Lorentz factors of the forward 
shock, reverse shock and the contact discontinuity are denoted by Tf{t), Tr(t) and Tc{t), 
respectively, and satisfy the relation F,. < Fc < Fj. 

The thin shell approximation assumes that the shocked layers are uniform. Then 71 = 
72 = Fc. For the situations envisaged here the forward shock can be considered ultra- 
relativistic, F/ >> 1. The jump conditions at the forward shock then yield F/ = v^Fc, 
and ^ 

Pi = gP^r^. (1) 

The reverse shock, on the other hand, cannot be considered ultra-relativistic in general 
and, therefore, a complete treatment is required. The jump conditions at the reverse shock, 
derived in appendix [Aj yield the relations 

p2 = Pehiq2,qe), M = beKq2,qe), P2 = Pe/(g2,ge,Cre), (2) 

subject to the condition 

^(g2,ge,c^e) = 0, (3) 

here q2 = (Fc/Fr)^, q^ = (•y^/Tr)'^, and the functions /, h and are defined in Eqs. flAlTD - 
flAlQp . Pressure balance at the contact, viz., pi = p2+b2/2, yields (4/3)pjF2 = p^f+blh'^/2 = 
Peif + c'"e^^/2), where Eqs. ([T]) and have been employed. Dividing the latter equation 
by (4/3)pj7g, noting that T'^/'j'^ = q2/qe, and defining 

G(r,t) = 3pe/(4p,,7e'), (4) 

one finally arrives at 

q2/qe = G{Rr,t)[f{q2,qe,(Te) + a e{Rr , t) h'^ {q2 , qe) / 2] . (5) 

The functions G{r, t) and <Je{r, t) in Eq. (|5]) are computed just upstream of the reverse shock, 
at r = Rrit), where 

Rr{t) =RrO + Vrdt' = RrO + j^^ (^1 " dt' (6) 

is the trajectory of the reverse shock, given to order 0(r^^) by approximating the velocity 
of the reverse shock as K — 1 — 1/2F^. Here t = to is the initial time of impact. Equations 
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([3]), ([5]), and determine the evolution of the variables q2{t), qe{t), and Rr{t) once G{r,t) 
and ae{r,t) are specified. 



The shock compression ratio is given by 

P272 (?2 + l)(ge-l) 



Pel'e (g2-l)(ge + l) 



(7) 



where 7g = 'je^ri^ — VeVr) and 72 = ■y2^r{^ — V2Vr) are the Lorentz factors of the unshocked 
and shocked ejecta, as measured in the shock frame. For an unmagnetized, relativistic shock 
Qe » 1, q2 — 2 and r = 3 as it should. The fast magnetosonic Mach number of the fiow 
just upstream the reverse shock is given by 

M,(t) < 



with Mg = (7g^ — ly^"^. The reverse shock exists as long as Ma > 1. For Ma < 1 the 
compression of the shell is communicated by a magnetosonic wave that propagates from the 
contact discontinuity backwards in the fiuid rest frame. 



2.1. Unmagnetized shell 

The limit of very low magnetization simplifies to de = in the above equations. To 
illustrate the properties of the solutions we examine two situations. The first one is that 
of a uniform shell, dr'^e = drPe = 0. For a non-expanding shell mass conservation implies 
Pe = Peo(t/to)~'^ and G{Rr,t) = (3peo/4aj)7~^(t/to)'^~^- For sufficiently small values of G the 
reverse shock is relativistic, yielding ~ 2r,^, 7^ >> F^, or equivalently q2 — 2, » 1. 
Eq. (lAlSp with cTg = 0, a = 4, ^2 = 2, qe » 1 gives f{q2, qe) — qel^- Substituting the above 
results into Eq. ([5]) one finally obtains 

Tc{t) =^l/\pJ^a,f'\tltot-'^'\ (9) 
recovering earlier results (Sari & Piran 1995). 

Our second example is the self similar ejecta considered by Nakamura and Shigeyama 
(2006; NS06), and discussed in appendix [Bl In this case 7e(x) = '^oX~^^'^ and the density 
profile is unrestricted. Here x = (1 + 2Fq)(1 — r/t) is the self-similar parameter of the 
unshocked shell and Fq is the initial Lorentz factor at the shell's head (x = 1) just before 
impact (at t = to)- The location of the reverse shock is given by 



Xr{t) = {l + 2Tl){l-Rr/t). 



(10) 
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Fig. 2. — The evolution of qe = {'je/Tr)^ (left panel) and q2 = (Tc/Tr)^ (right panel) for 
adiabatic index 7 = 5/3, k = (a constant density medium), n = 2, and different values of 
the parameter Gq = 3ae/4ajrg. 

Differentiating the latter equation, using dRr/dt = 1 — l/2r^ from Eq. (E]), and omitting 
terms of order Fq ^ and higher, gives 

^(tXr) = {To/Trr = (11) 
dt %{Xr,t} 

The initial conditions at t = to are Xr = 1, 7e(Xr) = To- For illustration we adopt pe = 
0'e{t/to)~^X"'^'^ for which 

G{xr,t) = {3aj4a.Tl){t/to)'"Yr^''/'. (12) 

Equations ( ITTl) . ([3]) and (j5]), with CTg = and G given by ( |T2l) . are solved simultaneously to 
yield Fc(t), Tr(t) and x,.(t). 

A particular solution can be sought for which q^, q2 are constants. Equation (fTTj) with 
7e(Xr)/ro = Xr readily yields 



Xrit) = {t/t,Y^^-'\ F, = q:'T,{t/t,Y'-^^^'\ F, = y^F,. (13) 
Eq. ([5]) with dqe/dt = dq2/dt = and Ue = implies dG/dt = 0, and using (fT2l) one finds 

9e - 1 = ^ir, 14 
n + 2 
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in accord with the self-similar solution derived by NS06. The value of q2 is computed 
numerically from Eq. ([3]). For any initial condition different than the value given in Eq. 
f lT^ self-similarity is broken and qe, q2 must evolve with time. However, it is found that 
for the range of values for which the NS06 solution is applicable the dynamics of the shell 
eventually approaches the self-similar limit, as naively expected since there is no scale in 
the problem. This is demonstrated in Fig [21 where solutions for qe and q2 are plotted 
for different initial conditions. As seen, the evolution quickly converges to the self-similar 
solution, whereby qe = 2.5 and q2 = 1.26 for the choice of parameters in this example {k = 0, 
n = 2). 



2.2. Magnetized Shell 

The structure of an impulsive high-o" shell has been computed recently by GKSIO and 
Lyutikov (2010a,b). They considered a situation in which a shell, initially at rest, is ex- 
panding into vacuum by pushing against a conducting wall. The shell is assumed to be 
uniform initially with a density po, magnetization o"o = ^o/Po >> 1? and finite width tq. 
For a spherical shell the total energy is E = (47r/3)rQpo<7o- The dynamics of the shell fol- 
lows several phases. At t = the shell starts accelerating, a simple rarefaction wave forms 
and propagates from the the shell's head towards the wall. At t = vq the rarefaction wave 
reaches the wall, reflects, and starts propagating back towards the head of the shell. At 
this point the shell loses causal contact with the wall. The major fraction of the energy 
(as measured in the Lab frame) is contained in a layer of thickness ~ 2ro bounded by the 
reflected wave (at the rear end) and the shell's head (see appendix IB.2p . The structure of this 
layer can be described by a self similar solution. The rear boundary of the shell accelerates 
as 7^^, ~ \/2(To^^(t/ro)^/^ (corresponding to the local magnetosonic speed in the fluid frame) 
and the fluid at the boundary as F^ ~ {aQt/2rQy^^. The head moves at a maximum Lorentz 
factor Fq = 2(Tq. In terms of the self-similar variable 

X = 8al[l-ir-2ro)/t], (15) 

the solution for the magnetization, density and Lorentz factor is given approximately by (see 
appendix [B]) 

^eix) = Pe{x)=Po^{t/r,)-\ 7e = 2aoX-^/^ (16) 

This solution is applicable above the front of the rarefaction wave, at % < x*(t), where 

xAt) = [1 + (SaoVoA)^/^/' ^ 8al{t/ro)-' (17) 
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Fig. 3. — Left panels: time profiles of the Lorentz factors of the unshocked ejecta 7e, con- 
tact discontinuity Tc and the reverse shock F^, for adiabatic index 7 = 5/3, k = and 
3po/(8o"Qaj) = 0.1 (upper window), 3po/(8c'"Qaj) = 10^ (lower window). Right panels: the 
evolution of the shock compression ratio r (see Eq. [7]), the ratio of kinetic-to- magnetic 
pressure /3 = 2^2/^2; the fast magnetosonic Mach number Ma of the fluid upstream the 
reverse shock (see Eq. [8]). 



(see Eq. |B21j ). The head of the shell is located at x = 1- 

To study the effect of the ambient medium on the evolution of the shell during the 
acceleration phase we assume that the structure of the unshocked shell is given by f|T6|) . We 
then obtain 



3po 



x'/Ve(Xr)(t/ro) 



fc-2 



The trajectory of the reverse shock is governed by the equation (see Eq. [TT]) 

^2 



N _ QeiXr)^ _ QeiXr) 2/3 



;i8) 



(19) 



At t = Xr = 1 and pe = o"e = 0. A strong shock forms initially at the shell's head and 
quickly propagates inwards (in the fluid rest frame). This brief initial phase is an artifact 
of out initial conditions that implicitly assume that impact with the external medium starts 
only after reflection of the rarefaction wave by the central source. As the density and 
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magnetization increases the shock weakens and the contact discontinuity accelerates until 
reaching a terminal value (Figs|3]and|l]). This phase lasts for a very short time. Subsequently, 
the contact discontinuity starts decelerating if /c < 2 or maintain a constant speed if = 2. 
For k > 2 the shell will eventually approach free expansion with — 2ao. For » 1, 
q2 » 1 we have from Eq. (]A17p h!^ ~ qe/(l2- Equation ([5]) with / = then admits the 
solution 

f 3poo-o\ ^^"^ 



rc(t) 



(t/ro 



l(fc-2)/4 



(20) 



V 8«i / 

which is viable as long as Fc < Fq = 2aQ. At early times, when q^ and q2 are of order unity 
the solution must be obtained numerically. Note that for k = 2 Eq. f l2U]) reduces to Eq. 
(47) in Lyutikov (2010a), that gives the Lorentz factor of a one- dimensional planar shell 
expanding in a constant density medium. This is expected since for the case considered here 
(a cold shell with a purely toroidal magnetic field) the flow equations in spherical geometry, 
(IBip - (lB4p . reduce to those in planar geometry upon the change of variables: p — )■ r^p, b — )■ rb. 
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Fig. 4. — Left panel: time profiles of Fc for 7 = 5/3, k = 0, and different values of Gq, where 
Co = 3po/(8o"Qaj). Right panel: the corresponding trajectories of the reverse shock Xr(t)- 
The red dotted lines mark the location of the rear boundary of the shell (the front of the 
reflected rarefaction wave) x*(^)) o"o = 10 (lower curve) and ctq = 30 (upper curve). The 
intersection of the two trajectories Xr(t) and x*(^) gives the time at which crossing of the 
accelerating shell by the reverse shock has completed. 
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Fig m delineates solutions for = (a constant density medium) and different values 
of Go = 3po/(8crQaj). For Gq » 1 the effect of the ambient medium is initially small; the 
reverse shock is confined at the head of the shell, viz., Xr ~ 1 (see right panel of FigH]), and 
Fc ~ Fq and declines slowly. For sufficiently high ctq the reverse shock will eventually start 
accelerating inwards in the shell frame before crossing of the shell has completed, so that 7e 
becomes significantly larger than F,. and ge >> 1, 92 >> 1- As seen in the left panel of Fig. 
m the evolution of the contact discontinuity then approaches Fc oc as expected from 

Eq. ( I20p . This transition occurs at a time t = tdec, where 

at which G{Rr, tdec) in Eq. (|5]) approaches unity. Here E = po<7o4vrrQ/3 is the total energy of 
the shell. The transition will occur before complete shock crossing, that is at 1 < xi^dec) < 
X*(W), if tdec < 8agro, or 



1 / 9E 



l/{8-2fc) 

(To > (Tcrt = ^ ( ) • (22) 

The latter scaling can be derived also from energy considerations. The energy accumulated 
inside the shocked ambient medium layer (i.e., between the forward shock and the contact 
discontinuity) at time t is approximately Ef = 47rFQajrQ(t/ro)^~'^/(3 — k). Deceleration of 
the shocked layer will commence before the shell reaches the coasting phase ii Ef c:^ E at 
time t < SaQVo. With Fq = 2cro the latter condition yields 



1 

an > - 
" 2 



(3 - fc)EV/(^-^'') , , 

^ ^ , (23) 



327ir^a 



in rough agreement with (|22|) . For a burst of duration T = 3OT30 sec and total energy 
E = 10^^i?52 ergs, expanding in an ambient medium of constant number density rii measured 
in c.g.s units the condition f l2^ reduces to 

/ E \^^^ 

ao > (Jert = 90 , (24) 



where the initial shell width has been taken to be tq = cT. Consequently, sub-critical shells 
(i.e., (Tq < (Tcrt) will not be significantly affected by the ambient medium before reaching the 
coasting phase. High-sigma shells (ctq >> Ccrt), on the other hand, will experience significant 
deceleration of the head well before the onset of the coasting phase, and it is anticipated 
that the Lorentz factor of the shell will be limited to Fmax — '^cTch- The dependence of the 
terminal Lorentz factor on the initial magnetization ctq) obtained from numerical integrations 
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Fig. 5. — The dependence of the terminal Lorentz factor of the accelerating shell Tmax, 
measured at time of complete shock crossing, on the initial magnetization ctq, for k = 0, 
7 = 5/3 and ^52/(71)71^) = 1 (solid line), Er^^/iTl^rii) = 100 (dashed line). 

of Eqs. ([3]), ([5]) and ( |T9l) . is exhibited in Fig. ([5]). As seen, for high-sigma shells the scaling 
r„,ax ^ 2acrt = 180(E52/T|oni)i/^ derived above using heuristic arguments, is in quite good 
agreement with the numerical result. Note that the critical magnetization depends on 
the initial energy density of the shell, E/tq. Thus, for a given power L = E/T = cE/vq we 
have acrt oc L^^^Tq ^ implying a less restrictive constraint on the Lorentz factor for smaller 
(sub) shells. For all cases studied it is found that a reverse shock always exists (i.e., Mj^ > 1) 
in the acceleration phase, however, except for the very early stages of the evolution the shock 
is weak and magnetically dominated, that is /3 = 2^2/^2 << 1- Emission from the shocked 
ejecta is not anticipated in high sigma shells, at least not in the ideal MHD case. 

In case of a stellar wind {k = 2) the shocked layer, after a brief rearrangement phase, 
maintains a constant speed (Fig. [6]), as expected from Eq. (120|) . The reverse shock in 
this case quickly weakens and eventually dies away {Ma becomes smaller than unity). The 
communication with the contact discontinuity then proceeds via a magnetosonic wave. For 
a mass loss rate M^, = M_5l0~^ Mq yr~^ and terminal velocity = lO^v^s cm s~^ the 
ambient density scales as rij = 10^^ M_^{r / Ro)~'^v~l cm~'^, with _Ro = lO^^-Ris cm being the 
radius at the wind's base. With this parametrization Eq. fl20l) yields for the Lorentz factor 
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of the contact 



10 . (25) 




Fig. 6. — Left panels: time profiles of the Lorentz factors of the contact discontinuity 
and the reverse shock F^, for adiabatic index 7 = 5/3, k = 2 and 3po/(8croai) = 10'^ 
(upper window), 3po/(8crQaj) = 1 (lower window). Right panels: the evolution of the shock 
compression ratio r, the ratio of kinetic-to-magnetic pressure (3 = 2p2/h\, and the Alfven 
Mach number Ma- 

As explained above, in a uniform circumburst medium a shell with initial magnetization 
(Jo < (Tcrt will start decelerating only after reaching the coasting phase, when conversion 
of magnetic energy into bulk kinetic energy has nearly completed. Once approaching the 
coasting phase the shell starts spreading and its structure is altered (GKSIO). To compute 
the evolution of the system in this region we adopt the self-similar solution derived in GKSIO 
(see also appendix |B]) for the unshocked, coasting ejecta. As noted in appendix [Bl this 
solution is not fully self-consistent in vacuum, as it implicitly assumes a confining agent 
at the head. However, is can be matched self-consistently to a shocked layer through the 
jump conditions at the reverse shock, and so may provide a reasonable description for the 
interacting ejecta after its spreading. With this choice 7e = FqX"^^^, G{x,t) is given by Eq. 
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Fig. 7. — The evolution of qe = {'je/Tr)^ (left panel) and q2 = (Tc/Tr)^ (right panel) for 
adiabatic index 7 = 5/3, k = 0, n = 2, 3ae/4ajrQ = 0.1, and different magnetizations. 

( !T2|) (see appendix |B|) . and 

<Je{x,t) = a,x-^-/'^'\t/t,,c)-\ (26) 

where ah denotes the magnetization at the head (x = 1) of the freely coasting shell at time 
of impact, t = tdec- The energy density of the unshocked shell can be expressed as 

= pe(l + ^e)ll = aer2(t/W)-3;^W2-l)(l + CT^), (27) 

with (Tg given by f l26p . so that for CTg << 1 the choice n = 2 corresponds roughly to a uniform 
energy distribution (since the shell is thin). Inspection of Eq. f l26p reveals that for n > —2 the 
magnetization of the shell decreases with increasing x- Consequently, it is anticipated that 
the solution will eventually converge to that of NS06 discussed in § §2.1[ This is confirmed 
in Fig. [71 In practice the shell has a finite width A. At the onset of the coasting phase the 
width of the shell is still ~ 2ro, corresponding to Ax ~ {2rQ/rc)2Tl ~ 1, using ~ ctq 
and Fq ~ ctq for the coasting radius and terminal Lorentz factor (GKSIO and appendix |B|) . 
At this point the width of the shell starts growing and at the deceleration radius r^ec it is 
expected to be A ~ TjAc, corresponding to Ax ~ r], where ?7 is a fraction of the ratio rdec/fc- 
Fig IS] exhibits solutions for a shell of width Ax = 10 interacting with a constant density 
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medium {k = 0). Although cTg < 1 is expected in the coasting phase, we present also a 
case with ah = 10 to elucidate the general behavior of the system. As seen, for cxh < I the 
compression ratio of the reverse shock increases considerably and the pressure downstream 
becomes kinetic dominated before complete crossing of the reverse shock, at \n{t/tdec) = 0.5. 
We find this behavior to be quite robust for self-similar ejecta with n > 0. 




Fig. 8. — Evolution of the Lorentz factor of contact discontinuity (left), ratio of kinetic-to- 
magnetic pressures /3 = 2^2/^2 (middle), and shock compression ration r (right), for a shell 
of initial width Ax = 10, k = 0, n = 2, and different values of the magnetization ah- The 
time at which reverse shock crossing is completed is ln(t/trfec) = 0.5 



2.3. Multiple shell model 

The rapid variability observed in many GRBs is commonly attributed to ejection of 
many sub-shells of small width that collide at relatively large radii. Such a multi-shell 
scenario can be envisaged also for magnetically dominated ejecta. In that case collisions of 
shells should occur only after full conversion of magnetic energy into bulk kinetic energy is 
accomplished, and before the onset of deceleration by the ambient medium. As explained 
above the restriction on the dynamics of a single shell imposed by the surrounding medium 
(Eq. [21]) is less stringent for a shell of small width, so that in principle it is anticipated that 



- 15 - 



for sufficiently rapid ejections the resultant shells may accelerate to the maximum Lorentz 
factor, r ~ (Jo ( at which the initial magnetic energy is fully converted), before a substantial 
fraction of their energy is dissipated at the forward shock. However, in order that the 
shells will collide only after reaching the coasting phase, the time interval between successive 
ejections of shells must be large enough, otherwise the shells will merge to form a steepening 
magnetosonic wave train with a characteristic size of the order of the engine's life time. 

To estimate the duty cycle required for shells to collide at the coasting phase, consider 
two shells of initial width /j and magnetization cTjo {i = 1, 2) expelled into vacuum. The 
shells will colhde provided (T20 > ctiq. Suppose that the ffist shell is ejected at time to 
and the second one at to + At. Now, the rarefaction wave of the ffist shell accelerates as 
7u;i = V^alo^{t/liY^^ (see Eq. |B19] ). Let 762(^2) = 2cr2oX2 ^''^ denotes the Lorentz factor 
of some point (that is, a fixed value of the self-similarity parameter X2) above the front of 
the rarefaction wave of the second shell. The relative velocity between that point and the 
rarefaction front of the first shell is 



4/3 / . \ 1/3 
'-'10 ' 

„2 



{2i 



1 

27^1 27,^2 4 
The two parts will catch up at time 

tcoll = / ^ Ahati' / -dy, (29) 

Jo Ve2-V^i Jo 1 - (2af/V7e2)l/'/' 

where yi = At/li. For 1a\^oV\'^ << 7e2 the result is approximately 

t,ou ^ ?>h{a^o^t/hY/\ (30) 

The coasting time (radius) of the ffist shell is t^ost — Sctiq^i- Thus, collision will occur after 
the ffist shell has reached the coasting phase provided tcou > icost or 

{At/h) > (8/3)3/V;/l (31) 

For a 10 » 1 Eq. fl3T]) imphes an unlikely small duty cycle. The above discussion suggests 
that within the framework of ideal MHD intermittent ejection of magnetically dominated 
outflow leads ultimately to one impulsive shell of size roughly equals the life time of the 
system. 



3. Conclusions 



We considered the interaction of a relativistic magnetized shell with an ambient medium, 
focusing on the case of an impulsive, high-sigma shell. We find that for values of the initial 
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magnetization ctq larger than the critical value given by Eq. fl22|) the evolution of the system 
is significantly altered by the ambient medium well before the shell reaches its coasting 
phase. For such high sigma shells a major fraction of the explosion energy is dissipated 
behind the forward shock by the time compression of the shell by the ambient medium is 
communicated to the accelerating rarefaction wave. The maximum Lorentz factor of the 
shell is then limited to values well below ctq (see Fig. [5]). Such episodes are expected to 
produce a smooth, relatively fast rising slowly decaying (power law) light curve, even in a 
multi-shell scenario. Events like GRB080916C and GRB090510 (Abdo et al. 2009a,b) are 
not easily accounted for by the impulsive high-sigma shell model. If extracted magnetically, 
such outfiows may require magnetic dissipation beyond the ideal MHD limit, as may occur 
in e.g., a striped wind model. Intermittent ejection doesn't seem to help, as unlikely small 
duty cycle is required in order for shells to collide after reaching the coasting phase. 

Lower sigma shells start decelerating only after reaching the coasting phase. The proper- 
ties of the reverse shock then depend on the structure of the coasting shell. For a self-similar 
shell with a reasonable density profile the magnetization decreases inwards and it is antic- 
ipated that the reverse shock will become strong before complete crossing. The evolution 
of the shocked layers is shown to quickly approach (for shells with a < 1) the self-similar 
solution derived by Nakamura and Shigeyama (2006). The slow acceleration of high sigma 
shells, (F oc r^/^), relative to a radiatively driven outfiow (F oc r), implies a smaller optical 
depth at the radius of collision. This can alleviate the need for extremely high Lorentz 
factors for most GRBs, to avoid strong absorption of the highly variable prompt emission. 

I thank Yoni Granot, Yuri Lyubarsky, Maxim Lyutikov, and Udi Nakar for useful com- 
ments, and the anonymous referee for a detailed report. This work was supported by an ISF 
grant for the Israeli Center for High Energy Astrophysics. 



where p is the proper density, p = p + pf, is the sum of kinetic pressure p and magnetic 
pressure ph = = 6^/2, h = e + p/p + b'^/p is the generalized specific enthalpy, is the 

fluid 4-velocity and g'^" is the metric tensor. The 4-vector of the magnetic field is defined in 
terms of the electromagnetic tensor F'^'^ as 



A. Jump conditions at the reverse shock 



The stress-energy tensor of a magnetized fluid takes the form 



T^^ = phu^'u^ - g^^p - b^W, 



(Al) 




(A2) 
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where rjap-ys is the Levi -Civita tensor. The jump conditions are obtained from integration 
of the flow equations, 



0, 



d^T'^'' = 0, 



d^ib^u" - b'u^) = 0, 
across the shock surface iIj{x^) = r — R{t, 6, (p) = 0, whereby we have 
[pu^jn^ = 0, [b^u'' - b^u^] = 0, [T'^^] = 0. 



(A3) 
(A4) 
(A5) 

(A6) 



The square brackets denote the difference of the enclosed quantity across the shock front, 
and 



(A7) 



Rr{t) 



is a 4-vector normal to the shock front. For the reverse shock with R{t, 6, (p 

n^ = r,,(-K, 1,0,0), 

here Vr = dRr/dt is the velocity of the reverse shock and = (1 — V^)^^^'^. Equations ( lAll) 
(IA6I) and (IA8I) yield for the quantities defined in Fig. [1] 

beleiVe - K-) = ^272(^^2 " K-), 

Pele{Ve - K) = P272(^^2 - K), 
PehellVe{Ve - K) + Pe = P2/i272^2(f2 - K) +P2, 
pJlelliVe - Vr) + PeK = P2h2ll{v2 - K) +^2^- 



(A9) 
(AlO) 
(All) 

(A12) 



We assume a cold ejecta, pe = 0, and adopt /i2 = P2 + 0.^2 where a = 7/(7 — 1) and 7 is 
the adiabatic index, for the shocked ejecta (see FigH]). Then, to order 0(F7^) the solution 
of Eqs. flA9l)- (lAT2|) is given by 

*(92,ge,0-e) = 0, 
P2 = Peh{<l2, Qe 
b2 = be{p2/Pe) = beh{q2,qe 
P2 = Pe/(g2,ge,Cre 

with 



(A13) 

(A14) 
(A15) 
(A16) 



h{q2,qe 



<l2 



fiq2,qe,(Te 



{qe - 1) 



a(g2 - 1) + 2 



1 - \/q2lqe + 



qe<ye _ <yeql{qe - 1) 

ge - 1 qe{q2 - 1)^ 

qe - 1 



(A17) 
(A18) 



^(^2, ge, CFe) = 2qe{l + (Te - ^Jq^j ^J(U) + (ct " 2)^20", 



52-1 
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+ (1 + a,)[a(g2 - 1) + 2] ( ^-^-^) - aa^q, (^-^ ] , (A19) 



'?2 + 1 / \qe 

where q2 = (Tc/Tr)'^, = {•ye/'^r)'^, and Ue = hi/ p^. The relations flA17P - flA19l) generahze 



those derived by Kennel & Coroniti (1984) for a high Alfven Mach number shock, corre- 
sponding to ge >> <l2 in our notation. In this limit Eqs. flAlSp and f lA19p with 7 = 4/3 
reduce to 

'-^Vf^V-^Uo. (A20) 



q2 \q2 - I J \l + CTe 

This condition can be expressed in terms of the 4- velocity of the downstream fluid, measured 
with respect to the shock frame, u'2 = rr72(f2 — K) = {12 — 1)/2-y/^ + 0(r~^), as 



1 + - Uu'i - = 0, (A21) 

which is equivalent to Eq. (4.10) in Kennel & Coroniti (1984). The solution for U2 is given 
by Eq. (4.11) in the same reference. 



B. Self-similar ejecta 

We consider a spherically symmetric, magnetized ejecta with a purely toroidal magnetic 
field, viz., 6^ = (0,0,0,6e)- The flow equations flASp - flASp then reduce to 

d{pele) I d 2 ^ n fT^T\ 

a(6e7e) 1 

r,. + -TrirbeleVe) = 0, (B2) 

ot r or 

|ln(p,/p2) =0. (B4) 

We seek self-similar solutions that are separable, to order 0(r^^, in the variables r = ln(t) 
and 

X = {1 + 2(m + l)r2}[l - (r - ra)/t], (B5) 

with the front of the expanding ejecta located at x = 1; and Fq = 7e(x = 1) = At'™-/"^ = 
^g-mT/2_ -yYe adopt the following parametrization of the fluid variables: 

7e = ^Igix). (B6) 

K = boe-P^B{x), (B7) 
p'e = Pele = p'oe-'^Hix). (B8) 
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Transforming from the coordinates (r, t) to (x, t) and using the relations 

tdt = dr + [(m + l){2Tl -x) + (B9) 

tdr = -[I + 2{m + l)Tl]d^, (BIO) 

(Bll) 

one obtains, upon substitution of Eqs. (1B6P -( 1B8|1 into the flow equations (IBip - (lB4p . 

= m(l - gx) - 4ae(l - p - m/2), (B12) 

= ^ _ (1 _ _ _ 2ae(l - p - m/2), (B13) 

g dx 

7^ = m - (2 - g)(l - ^7x) + (Y^^K2 + p - 1 + (2 - g)<7x]. (B14) 
to order 0(Tq'^), where A = {m + 1)[(1 — gxY ~ ^(^edxl- 



B.l. Freely expanding ejecta 

For a freely expanding ejecta d'je/dt = —m/2 + (m + l)(l/5' — x)9x\f9 = 0. The 
boundary condition g{x = 1) = 1 implies m = and gix) = X^^- Equations flB12p - flB14p 



give p = 2, q = 3, B(x) = X The function H{x) is unrestricted. To summarize, the 
solution can be expressed as 



7, 



'-Tlx-\ K = h,e-^^x-"^ p', = a,V,e-^^H{x). a, = aoe-^x'^^'H-^x), (B15) 



e 



with bo denoting the value of 6e at r = x = 1- The discussion in NS06 suggests that a 
power law density profile, pe oc 7~", corresponding to H = x^"'~^^^'^i provides a reasonable 
description of realistic ejecta. Note that magnetic and kinetic energies have different scaling, 
Pe7e ^ aiid 6g7g oc t~*, so that the enthalpy is not self-similar. Note also that magnetic 
energy is not conserved, viz., b^t^ oc t~^, implying a loss of Foynting energy from the front. 
This is a consequence of an implicit boundary condition at the vacuum-shell interface that 
assumes a confining agent. The loss of magnetic energy is then associated with a pdV 
work. For the interacting shell this solution can be matched with the shocked shell layer, as 
discussed in ^2.2[ 
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B.2. An accelerating high-a shell 

The solution describing an impulsive accelerating shell (GKSIO, Lyutikov 2010a) corre- 
sponds to the choice m = 0, p = 1, q = 2. Eqs. flB12|) -f lBT4|) then readily yield 



-e(x,r) = ^^^^, (B16) 

independent of r, and 

9{X) = X-^'\ (B17) 

B{x) = Hl^g = a,[x)- (B18) 

The requirement 7e = VqxT^^^ > 1 formally implies X < Xo = ^o^^ and Fq < 4cro- However, 
the above solution is applicable only well above the magnetosonic point, defined by the 
condition (Tg = ■y'^v'^ ~ 7g. A full treatment (GKSIO) gives Tq = 2(Tq where (Tq is the initial 
magnetization of the shell. To the order at which we are working the self-similar variable is 
then given by (fT5|) and the solution by (fT6|) with po denoting the initial density of the shell. 

We choose = = 2ro in (1B5I) . where vq is the initial shell width and is the 
magnetosonic radius. Then at r = (corresponding to t = r^) the front of the refiected 
rarefaction wave is located at x = ScTq. The local velocity of the front, as measured in the fiuid 
frame, is v'^ = \/ a^j (1 + cTg), and the corresponding Lorentz factor is 7^ = (1 — f^)~^/^. The 
Lorentz factor of the wave front in the lab frame is obtained upon a Lorentz transformation: 

1^ = lel'Jl + VeV'J ~ 27ev/^. (B19) 

The trajectory of the wave front is governed by the equation dr^/dt = — I — 1/27^, which 
can be translated into 

J 2/3 

Sta.) = (8../7.)^ = ^.XP. (B20) 

equality h( 
equation reads 

X.{t) = [1 + {Salro/tf'f/^ (B21) 

At t « Sa^ro we have x* ^ 8al{ro/t) and 7e(x*) ^ (ao^/2ro)^/^ 7«; - V2al^\t/roy/'^ in 
agreement with GKSIO and Lyutikov (2010a). 

The energy density is 

= pe(l + a^H - ^(1 - X-'^nt/ror. (B22) 



1/3 

where the last equality holds at x* >> 1 for which 4o"e — x* • The solution of the latter 
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Using r'^dr = {t^ /8crQ)dx we obtain for the total energy contained between the reflected 
rarefaction wave and the head 



AE = / T^Hnr^dr = npocrorl{t/8al) / (1 - x^'^^dx ^ vrpo^or|| = {3/A)E (B23) 



J Ji 
independent of time, where E is the initial energy (i.e., the explosion energy). 
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